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We consider the simulation of the dynamics of one nonlocal Hamiltonian by another, allowing 
arbitrary local resources but no entanglement nor classical communication. We characterize notions 
of simulation, and proceed to focus on deterministic simulation involving one copy of the system. 
More specifically, two otherwise isolated systems A and B interact by a nonlocal Hamiltonian H 7^ 
Ha + Hb. We consider the achievable space of Hamiltonians H' such that the evolution e^'^ ' can 
be simulated by the interaction H interspersed with local operations. For any dimensions of A and 
B, and any nonlocal Hamiltonians H and H' , there exists a scale factor s such that for all times 
t the evolution e~*^ "* can be simulated by H acting for time t interspersed with local operations. 
For 2-qubit Hamiltonians H and H' , we calculate the optimal s and give protocols achieving it. The 
optimal protocols do not require local ancillas, and can be understood geometrically in terms of a 
polyhedron defined by a partial order on the set of 2-qubit Hamiltonians. 



I. INTRODUCTION 

A. Motivation 

Like the mythical lovers Thisbe and Pyramus, Alice and Bob wish to be forever in each other's company, a situation 
described physically by some many-atom interaction Hamiltonian H' . Unfortunately their parents disapprove, and 
have built a massive wall to keep the youngsters apart. Fortunately there is a small hole in the wall, just big enough 
for one atom of Alice to interact with one atom of Bob via the two-atom interaction Hamiltonian H (Fig. 1). Can they 
use this limited interaction, together with local operations on each side of the wall, to simulate the desired interaction 
H'l Yes, if they are patient, because any nontrivial bipartite interaction can be used both to generate entanglement 
and to perform classical communication. Therefore they can use H , along with local ancillary degrees of freedom 
on each side of the wall, to generate enough entanglement, and perform enough classical communication to teleport 
Alice's entire original state to Bob's side. Now that they are (virtually) together, Alice and Bob can interact to their 
heart's content. When it is time for Alice to go home, they teleport her back to her side, in whatever entangled 
state they have gotten themselves into, again using H to generate the needed entanglement and perform the needed 
classical communication. So, by the time they get to be old lovers, Alice and Bob can experience exactly what it 
would have been like to be young lovers, if they are still foolish enough to want that. 
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FIG. 1. Thisbe (L) and Pyramus, separated by a wall, through which they can only interact by a 2-atom Hamiltonian H. 
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A more practical motivation for studying the ability of nonlocal Hamiltonians to simulate one another comes from 
quantum control theory [y , in particular the problem of using an experimentally available interaction, together with 
local operations, to simulate the evolution that would have occurred under some other Hamiltonian not directly 
accessible to experiment. A more mathematical motivation comes from the desire to parameterize the nonlocal 
properties of interaction Hamiltonians, so as to characterize the efficiency with which they can be used to simulate one 
another, and perform other tasks like generating entanglement Ig,^ or performing quantum computation y-[7|. This 
parallels the many recent efforts to parameterize the nonlocal properties of quantum states, so as to understand when, 
and with what efficiency, one quantum state can be converted to another by local operations, or local operations and 
classical communication. It is not difficult to see, by the Pyramus and Thisbe argument, that all nonlocal Hamiltonians 
are qualitatively equivalent, in the sense that for any positive t' and e, there is a time t such that t' seconds of evolution 
under H' can be simulated, with fidelity at least 1— e, by t seconds of evolution under H, interspersed with local 
operations; but much work remains to be done on the quantitative efficiency of such simulations. 

In this paper we derive bounds on the time efficiency with which one Hamiltonian can simulate another using local 
resources. In the case of two interacting qubits, we show that these bounds are optimal. The structure of the paper 
is as follows. In Sec. ||, we define the allowed resources and the type of simulation we consider. In Sec. |l|, we prove 
some general results on the type of simulation we consider along with some examples. In Sec. P^, we define our goal 
and summarize our main results for two-qubit Hamiltonians, that are proved in Sees. |V| and]Vl|. Some discussions 



and conclusions, and more auxiliary results can be found in Sec. VII, VIII and Apps. |A[]D| . We first describe in more 
detail some related results. 

B. Related work 

The qualitative equivalence of nonlocal Hamiltonians noted above, and the use of interaction as an infinitesimal 
generator of entanglement, was already noted several years ago S. These discussions also considered the question 
of interconverting discrete nonlocal primitives, such as nonlocal gates, shared EPR pairs, and uses of a classical bit 
channel. More generally and quantitatively one may ask, given a nonlocal Hamiltonian Hj^g ^ Ha + Hg , what is the 
optimal efficiency with which it can be used, in conjunction with local operations, 

• to generate entanglement between A and B 

• to transmit classical or quantum information from A to B, or vice versa 

• to simulate the operation of another nonlocal Hamiltonian H' . 

A partial answer to the first question, for two-qubit Hamiltonians, was given by Ref. B. The current work is a 
continuation of previous efforts to study the efficiency simulating one Hamiltonian by another. 

Hamiltonian simulation has been considered in the context of quantum computation @- 7] j|9|-[ll|] . In these works the 
system consists of n qubits, with some given pairwise interaction Hamilton ian. I n Refs. \i-^, the given Hamiltonian 



was a sum of az ® ctz interaction terms between distinct qubits (see Sec. |III C for definitions) and the goal was to 



simulate a particular one of these terms. This was extended in Refs. [[7 yiO|Jll[| to arbitrary pairwise interactions, in 
both the simulating and the simulated Hamiltonians. In these papers the main concern was to obtain methods for 
simulation, and therefore upper bounds on the resources as a function of n. 

Independent results on optinfizing the time used of a given Hamiltonian for performing certain tasks are reported in 
Refs. [p p^Jl^ . Reference B gives a necessary condition for simulating one n-qubit pairwise interaction Hamiltonian 
by another, and gives a necessary and sufficient condition for simulation with a particular given Hamiltonian. Time 
resources for simulating the inverse of a Hamiltonian are discussed in Refs. |^,|l^,|l2[. Reference ||l^ considers simulat- 
ing a unitary gate using a given Hamiltonian and a set of controllable gates in the shortest time. A general framework 
is set up in terms of Riemannian geometry. A time optimal protocol is obtained for the specific Hamiltonian ct^ ® u^ 
in the 2-qubit case. 

Finally, some more recent results have appeared since the original posting of this preprint, extending it and related 
work in various ways @-0 [Ill-ill- 

II. SIMULATION FRAMEWORK 

In this section we describe our framework of Hanfiltonian simulation, i.e. the rules under which the simulation is 
to be performed. We also describe other possible frameworks and their relations to the one we adopt. 



A. Available resources 

Let H and H' each be a nonlocal Hamiltonian acting on two isolated systems A and B. We consider the problem 
of simulating H' by H using unlimited local resources. These include instantaneous local operations and uncorrelated 
local ancillas of any finite dimensions. It is also necessary to allow some initial classical correlation - Alice and Bob are 
assumed to have agreed beforehand on their time and spatial coordinates and the simulation protocol to be followed. 
Besides this, no other nonlocal resources are allowed, neither prior entanglement nor any form of communication 
beyond what can be achieved through the interaction H itself. Our goal is to minimize the time required of the given 
Hamiltonian H to simulate another Hamiltonian H' . This will be defined more formally in Sec. [rv| . 

Note that either the simulating or the simulated system or both can be given the freedom of bringing in local 
degrees of freedom (ancillas) and allowing interaction between each ancilla with the corresponding local system. 
Ancillas on the simulated system can make the simulated more powerful and therefore harder to simulate. Ancillas 
on the simulating system potentially make the simulation easier. We will allow ancillas on the simulating system. 



though they may not always help (Section VI) 



B. One-shot and deterministic simulations 

In this paper we only concern ourselves with protocols that are one-shot — i.e. operate on a single copy each of the 
simulated and simulating systems — and that are required to succeed with probability 1. 

More generally, a simulation can be "blockwise", in which iJ^" is used for the simulation of iJ"^", or in which 
H is time-shared among many copies of the system and the amortized cost is considered. A simulation can also be 
stochastic and fail with finite probability, in which case the expected cost is considered. 

C. Gate versus dynamics simulations 

One possible notion of simulation is that, given H' and t\ we simulate the final unitary evolution e~'^ * by 
composing local operations with elements in the one-parameter family {e~*^*}tJj The final evolution needs to be 
correct, but the intermediate evolution need not correspond to e~*^ * for < t" < t' . The efficiency, given by the 
ratio t/t' can depend on t' . For example, H can be used to generate entanglement and classical communication to 
bring A and B together by t deportation, apply e~'^ * , and teleport A back. Viewing the cost i as a function of 
t' , t does not increases indefinitely with t' , rather, it can be made constant after it reaches a sufficiently large value. 
As another example, if the nonlocal Hamiltonian H' — az 'S) ctz acts for time t' = 7r/2, the result is the unitary gate 
iuz ® (Jz, which is local, and requires no nonlocal interaction time at all to simulate. This type of simulation, with 
very different primitives, is much studied in the context of universality of quantum gates (composing a small set of 
available gates to obtain any desired unitary gate). More recently, simulation of a unitary gate using a fixed given 
Hamiltonian for a minimal amount of time and local manipulations was studied in Ref. |l3] and some partial results 
were obtained. From now on, we call this type of simulation "gate simulation" or "finite time simulation" . 

A natural direction to strengthen the above notion of Hamiltonian simulation is to require not only the end result, 
but also the intervening dynamics of if to be simulated. Intuitively, one might expect this to mean that the application 
of H, interspersed with instantaneous local operations, produces a trajectory that remains continuously close to the 
trajectory e~*^ * which one wishes to simulate. However, this is impossible in general, because the needed local 
operations cause the simulating trajectory to be discontinuous, agreeing only intermittently with the trajectory one 
wishes to simulate. Accordingly we adopt the following definition of dynamics simulation: The Hamiltonian H 
simulates the dynamics of H' with efficiency /i if Vt'>o, V(:>o the unitary operation e~*^ * can be simulated with 
fidelity > 1 — e by some protocol using H for a total time t' / ji and local operations. While this characterization 
may appear to have given up the idea of approximating the simulated system at intermediate times, in fact it has 
not, because it can be shown to imply the existence of a /i-cfficient "stroboscopic" simulation, which approximates 
the simulated trajectory arbitrarily closely not only at the begining and end, but also at an arbitrary large set of 
intermediate times. We discuss this and other simulation notions in Appendix H. We also show that the existence 
of a protocol for dynamics simulation is equivalent to the existence of one for simulating an infinitesimal time (see 



^ The evolution due to a Hamiltonian H is given by e '^'. Note the — minus in the exponent. 



Sec. [II A) which in turns imphes the abihty to create protocols for arbitrary finite times by appropriately rescaling 
and repeating the infinitesimal-time protocol, (see Appendix. |D| ). 

III. GENERAL RESULTS AND EXAMPLES 

Having defined the simulation framework, we derive some important general results and provide some examples of 
dynamics simulation, which motivate our main results and simplify some of the later discussions. 

A. Infinitesimal and time independent simulation 

First of all we show that dynamics simulation is equivalent to "infinitesimal simulation" , the problem of simulating 
the evolution of H' for an infinitesimal amount of time t' . On one hand, any protocol for dynamics simulation 
simulates the initial evolution, therefore is a protocol for infinitesimal simulation. On the other hand, iterating an 
infinitesimal simulation results in dynamics simulation. We restrict our attention to infinitesimal simulation from now 
on, and focus on the lowest order effects in t' . Note that this property may not hold for other types of simulation 
described in Appendix |A|. 

Infinitesimal simulation has a very special structure - the optimal simulation protocol is independent of the in- 
finitesimal value of t' . The proof is included in Appendix y. 

B. Local Hamiltonians are irrelevant 

A general bipartite Hamiltonian K can be written as, 

K ^KA®I + I®KB+^M^Jr|^®■qJ, (1) 

i'j 

where / denotes the identity throughout the paper, Ka, Kb are local Hamiltonians acting on A, B respectively, and 
{?7i} is a basis for traceless hermitian operators acting on each of A and B. We can "dispose" of the local Hamiltonians 
Ka and Kb by undoing them with local unitaries on A and B: 

In other words, K can be made to simulate its own nonlocal component. 

Likewise, any Hamiltonian can simulate itself with additional local terms. Therefore, given unlimited local resources, 
the problem of simulating an arbitrary Hamiltonian H' by another arbitrary one H reduces to the case when both 
are purely nonlocal. 

C. Possible inefHciencies in simulation 

Consider the simplest case of two-qubit systems. We introduce the Pauli matrices 



^ 1 ; ' ^^« ^ I i ; ' "" ~ V -1 ' ' ^^^ 

and the useful identity 

f/e^C/t = e^^^^' (4) 

where M is any bounded square matrix and U is any unitary matrix of the same dimension. 

As an example, let H = (7x®(Jx and H' = \{ax®o-x + o-y®(yy + <yz®<^z)- To simulate H' by H, let Ui = ^(cTaj+CTj,) 
and U2 ~ -7^{<^x + <7z), so that ay = UiaxUl and cr^ — U2crxUl. Using Eq. (H), it is easily verified that 

€-'"'*' = ( e-*^*'/3 ) X ( C/i (g) C/i e^'"''/^ Ul C// ) x ( C/2 (g C/2 e~'"''/^ Ul ®Ul) . (5) 



Conversely, we can simulate H with H': 

Note that the simulating of H' for a duration of t' requires applying H for a duration of t' whereas simulating H for 
a duration t requires applying H' for a duration of 3t. As the time required of the given Hamiltonian is a resource to 
be minimized, we see that some simulations are less efficient than the others. In this paper, we are concerned with 
the inefhciencies of simulation intrinsic to the Hamiltonians H and H' that arc not caused by a bad protocol. For 
example, we will show later that the inefficiency in the above example is intrinsic. 

D. Simulating the zero Hamiltonian — stopping the evolution 

In some applications, the given Hamiltonian H cannot be switched on and off. Simulating the zero Hamiltonian 
can be viewed as a means for switching off the Hamiltonian H ^-g] . This can always be done for any dimensions of 
A and B. 

First, let A and B be 2"-dimensional, and 

i/ = ^CijPi®Pj, (7) 

ij 

where i is a binary vector (Ji,?2, ■ • • ,*2n) that labels the n-qubit Pauli matrix Pi — a^^a"^^ ® ■ ■ ■ ® cr^^'^V*^". It is 
easily verified that 

^^ V PiMPi ^ trM — . (8) 

22n / ^ ' ' 2" 

i 

A protocol for simulating by iJ is given by, 

n,,{P®P,)e-^"'/^'-{P}®Pl) =. e-^Y.;iP^^P^)"iPl^Pl)+^i^') « ^-.*tr^/2-^ (9) 

in which the net evolution is just an overall phase to the lowest order in t. 

When A and B are d-dimcnsional, one can embed each of A and i? in a larger, 2"-dimensional system for n = [log2 d] 
to perform the simulation. Physically, this can be done on each of A and B, by attaching a qubit ancilla, extending 
the Hilbert space to 2d-dimensions, and applying the simulation to a 2"-dimensional subspace, such as one spanned 
by |z) (E)\0) for i = 1, • • • , d and \i) (g) |1) for i = 1 ■ ■ ■ 2"— d. Such simulation can also be done without ancillary degrees 
of freedom, and an alternative method based on Ref. [g3| is given in Appendix ^. 

E. Arbitrary but inefficient simulations 

We now show that any nonlocal bipartite Hamiltonian can be used to simulate any other, albeit with inefficiencies. 
In other words, for any H and H', operating H for time t can simulate the evolution of H' for time t' with t' /t > 0. 
This holds for any dimensions. We keep all definitions from the previous example in the following protocol. 

First, let A and B be 2^-dimensional, H = J2i\ cy-Pi (8) Pj and H' = J2i\ cJj-Pi ® Pj- Without loss of generality the 
coefficient for P^ (g) Pk is positive, i.e. Ckk > 0, where k = (0, 1, 0, • • • , 0) and Pk = CTz (X" / • • • /. It is known that 
for any Pi and Pj , there exist unitary operations C/yi in the Clifford group Q , such that 

C^iJ±^iC^ij± = ±^j (10) 

In other words, one can always transform any Pi to any other or to its negation. In our protocol, H simulates H' in 
two steps. First, H simulates Pk (Xi Pk by: 

Alice and Bob independently apply an averaging over all Pauli operators commuting with Pk, removing all operators 



except for / = Pq and Pk in each of their systems. The local terms can be ignored, following Sec. IIIB, Second, 
-Pk <8) Pk simulates H' by: 



where sg(a:) — x/\x\ if a; 7^ and wc omit terms with c^ ■ = . 

When A and B are d-dimcnsional, the simulation of sH' by H can again be performed in a larger 2" x 2" system. 
This method implies a lower bound on the maximum possible value of s, s > a r ioga ^1 '^^'■' ''^'^ . It is also possible to 

perform the simulation without ancillas. The proof is given in Appendix kV Other methods for such simulation were 
independently reported in |lq-|20||. 

F. Equivalent classes of local manipulations 

Under our simulation framework, Alice and Bob are given unlimited local resources. In this subsection, we show that 
they only need a relatively small class of manipulations. To facilitate the discussion, we introduce classes of operations 
C, that can be LU, LO, LU+anc, and LO+anc, to be defined as follows. LU is the class of all local unitaries that act 
on A (g) i?. LU+anc is similar, but acts on [A ® A') ® {B ^ B') where A' and B' are uncorrelated ancillary systems of 
any finite dimension. LO and LO+anc are similarly defined, with the unitaries replaced by general trace-preserving 
quantum operations. Note that the largest class LO+anc corresponds to what is most generally allowed under our 
simulation framework. 

We now show that LU+anc, LO, and LO+anc are equivalent under our framework. First, we show that LU+anc is 
at least as powerful as LO+anc. Any trace preserving quantum operation can be implemented by performing a unitary 
operation on a larger Hilbert space, followed by discarding the extra degrees of freedom (see, for example, Ref. p5[). 
The exact difference between LO+anc and LU+anc is that measurements and tracing are disallowed in the latter. 
However, these are not needed when simulating Hamiltonian in LU+anc, due to the following facts. (1) Measurements 
can be delayed until the end of the protocol, as operations conditioned on intermediate measurement results can be 
implemented unitarily. (2) In Hamiltonian simulation, the ancillary systems A' B' have to be disentangled from AB at 
the end of the simulation. Thus no actual measurement or discard is needed. These facts allow any LO+anc protocol 
to be reexpressed as an LU+anc protocol with pure product state ancillas, meaning that LO and LO+anc are no more 
powerful than LU+anc. Conversely, due to fact (2) above, any LU+anc protocol can be viewed as an LO protocol. 
Thus, we establish the equivalence between LO, LU+anc, and LO+anc. From now on, we focus on LU+anc protocols 
for full generality, and on LU protocols as a possible restriction. 

IV. FORMAL STATEMENT OF THE PROBLEM AND SUMMARY OF RESULTS 

Let H, H', A, B, A', B' be defined as before. 

Definition: H' can be efficiently simulated by H, 

H' <c H, (13) 

if the evolution according to e~*^ * for any time t' can be simulated by using the Hamiltonian H for the 
same time t' and using manipulations in the class C . 

Definition: H' and H are equivalent under the class C, 

H' =c H, (14) 

if H' <c H and H' <c H. 



Throughout the paper, we only consider LU +anc protocols following Sec. IIIF. We also restrict attention to H 



and H' that are purely nonlocal, following Sec. HIE 



An LU+anc protocol simulates H' with H by interspersing the evolution of H with local unitaries on A A' and BB' . 
More specifically, the most general protocol for simulating H' using H for a total time t is to attach the ancillas A' B' 
in the state \0a') ® \^b'), apply some C/i (g) Vi, evolve AB according to H for some time ii, apply C/2 ® ^2j further 
evolve AB according to H for time ^2, and iterate "apply Ui ® Vi and evolve with H for time ti" some n times. At 



the end, it applies a final Uf®Vf. The ii > are constraineqj by X]r=i ^J ~ *• Suppose the protocol indeed simulates 
an evolution for time t' according to H' . Then we can write 

= ( e-^^'*'|7/;) ) ® ( WA'B'itu • ■ ■ , tn)\QA') ® \Qb') ) , (15) 



where we have redefined Ui=i^2,---,n and 14=i, 2, •••,«, and |-0) denotes the initial state in AB. In Eq. (|l5|), e~*^*' acts 
on AB and implicitly means e~'^*' ® Ia'B' ■ The operator WA'B'{ti, • • • , in) describes the residual transformation of 
A'B', and can be chosen to be unitary since the operation on the left hand side of Eq. (M) is unitary. The problem 
we are concerned with can be stated in two equivalent ways: 

Optimal and efficient simulation: Let H be arbitrary. The optimal simulation problem is to, for each 
H' , find a solution {Ui}, {Vi}, {ti} of Eq. ( |l5[ ) such that t' /t is maximal. The efficient simulation problem 
is to characterize every H' which admits a solution for Eq. dl3) with t' = t, i.e. H' <LU+anc H. 

Definition: The optimal simulation factor sh'\h under class C of operations is the maximal s > such 
that sH' <c H. 

The optimal and efficient simulation problems are equivalent because inefficient simulation is always possible (see 



Section III). The efficient simulation problem can be solved by finding the optimal solution for each iJ' and charac- 
terizing those with t' /t > 1. The optimal simulation problem can be solved by finding the maximum s for which sH' 
is efficiently simulated. With this in mind, we may talk of solving either problem throughout the paper. 

We now summarize our results. We show in Appendix ^ that, in the infinitesimal regime, the most general 



simulation protocol Eq. (15) using LU+anc is equivalent to 



sH' = {Oa'\ ^ {Ob'\J2 P^ ^» ^ ^^ (H^Ia'B') Uj f/ 10^0 ® |0b') . (16) 

i 

In the LU case (without ancillas), Eq. (llW) reads 

sH' = ^p^U^®V,HU}® V} , (17) 

i 

where i = ii + • • • + t„, p^ = t^jt^ and s = j. Thus, the set {H' <lu -ff} is precisely the convex hull of the set 
\\J ®V H U'^ ^V'^^ when U and V range over all unitary matrices on A and B respectively. The linear dependence 
of ^-W on H is manifest in both Eq. (|l|) and Eq. (|l^. 

Our main results apply to the simulation of two-qubit Hamiltonians, and are summarized as follows: 

Result 1: Any simulation protocol using LU+anc can be replaced by one using LU with the same 



simulation factor. This will be proved in Section Vl. Thus, the four partial orders <lUi ^LU+ano ^lOj 
<LO+anc are equivalent for two-qubit Hamiltonians. 

Result 2: We present the necessary and sufficient conditions for B' <lu H, for arbitrary two-qubit 
Hamiltonians H and H' , and find the optimal simulation factor Sf{'\H and the optimal simulation strategy 
in terms oi {Ui} , {Vi} , {ti} . This will be discussed in Section 0. 

These results naturally endow the set of two-qubit Hamiltonians with a partial order <c- This induces for each 
H, a set {H' : B' <c B} which is convex: if B' <c B and B" <c B, pB' + (1 - p)B" <c B for any < p < 1. 
Our method relies on the convexity of the set {B' : B' <c B}, which has a simple geometric description, and in 
turns allows the partial order <c to be succinctly characterized by a majorization-like relation. The geometric and 
majorization interpretations offer two different methods to obtain, in practice, the optimal protocol and the simulation 
factor. 



Without loss of generality, a protocol with 5^"=! ti < t can be turned to one with ^^"=1 ii = i by simulating the zero 



Hamiltonian as described in Section 



[II 



V. OPTIMAL LU SIMULATION OF TWO-QUBIT HAMILTONIANS 

We will prove that <lu is equivalent to <LU+anc in the next section. In this section, we focus on LU simulations. 
We first adapt a result from Ref. [Q to reduce the problem to a smaller set of two-qubit Hamiltonians H and H' . 
Then, for any H, we identify the set {H' : H' <c H} with a simple polyhedron and obtain simple geometric and 
algebraic characterizations of it. The optimal solution for each pair of H and H' is derived. Finally, the problem is 
rephrased in the language of majorization. 

A. Normal form for two-qubit Hamiltonians 

The most general purely nonlocal two-qubit Hamiltonian K can be written as, 

ii: = ^ Afzj cr» CTj , (18) 

where the summation is over Pauli matrices i, j — x,y, z or 1,2,3 throughout the discussion for two-qubit Hamiltoni- 
ans. Let 

H = ^h,a^®ai, (19) 

i 

where /ii > /12 > |/i3| are the singular values of the 3x3 matrix M with entries Af^, and /13 — sg(detM)|ft,3|. We 
say H is the normal form of K. 

Theorem: Let H be the normal form of K . Then H =lu K. 

Proof: If the local unitaries U^ ® V^ and U ®V are applied before and after e~*^*, the resulting evolution 
is given by 



with 



K' = {U (i()V)K {U^ ®V^) 

ij 

= E ^^»^- ( E ^^'^' ) ® ( E ^^fc'^'^ ) (21) 

ij I k 

= ^ ( R^MS )^^al(g>ak ^ Yl ^ik '^i ® ^^ (22) 

Ik Ik 

In Eq. (Ell), R, S <E SO (3) since conjugating r • it by SU(2) matrices corresponds to rotating r by a matrix 
in S0(3j(and vice versa). Equation ( p2|) implies that K' — U ® VKW (g) V^ for some unitary U, V if and 
only if A/' = R!^ MS. In particular, there is a choice of R and S that makes K' = H: 

/100\ /ioo\ 

i?^=01 xOf, 5 = Ojx01 (23) 

\ det Oi y \ det O2 / 

where M = O1DO2 is the singular value decomposition of M, with Oi,02 G 0(3) and D = 
diag(ft,i, ft,2, 1^3 1). Thus K and H are related by a conjugation by local unitaries, which implies K =lu H- 

As suggested by the above proof, we define a few useful notations. 

Definitions We call the 3x3 real matrix Mij the "Pauli representation" of K, when M and K are related 
by Eq. (|l8|). We use Dk to denote a diagonal Pauli representation of K. 



Since any 2-qubit Hamiltonian is equivalent to its normal form, we assume iJ', H are in normal forms from now on. 
We now turn to LU simulation of H' by H. 



B. General LU simulation of normal form two-qubit Hamiltonians 



Recall from Eq. {njh in Section IV that the most general simulation using LU is given by 
sH' =piUx®ViH Ul ®Vl + ■ ■ ■ +PnUn®Vn H Ui® V^ 



(24) 



where s = t' jt. Following the discussion in Section V A , we only need to consider H = J2i^i '^i ® ^i ^'^'^ ^' — 
Si ^'i ^i ® '^i that are in their normal forms. The Pauli representation of {U ® V)H{W ®V'^) is given by RDhS for 
some R,Se SO (3). We can reexpress Eq. (p4) as 



sDh' = Pi RiDhSi 



Pn RnDnSn 



and h'^ > h'2 > 



(25) 
h'J. Without loss of 



where Ri, Si G S0{3). Since H and H' are in their normal form, hi > h2 > {h^ 

generality, we can make two assumptions. First, we can assume ^13 > 0: If ft.3 < 0, we can right-multiply Eq. (|2£ 

^ = diag(l,l,-l): 



by 



sDh'S = pi Ri{DhS)(SSiS) 



Pn Rn{DHS){SSnS) 



(26) 



in which SSiS £ S0(3), and DhS — diag(/ii, /i2, l^sl) is of the desired form. Thus, we can assume /13 > 0. Second, 
note that s^iii^ = o,Sfj'\aH = ^■5aff'|ff ■ The protocol is unchanged when Eq. ( |25| ) is divided by trDn — hi + h2 + /13. 
Therefore, without loss of generality, the normalization /ii + /12 + /13 = 1 can be assumed. 

Equations (|2J) and ( psj ) have a simple physical interpretation: the protocol partitions the allowed usage of i? (Dh) 



rt 



Vj {RkDSk), resulting in an "average Hamiltonian" H' {Dh> 



which is a convex 



into different Uk ® Vk H Ul 

combination of the Uk®Vk HUl® V^ (RkDSk). 

The Hamiltonians, represented by Dh' , that can be efficiently simulated (s = 1) correspond to the diagonal elements 
of the convex hull of {RDhS : R,S £ S0(3)}. We call this diagonal subset, which is also convex, Ch- Note that 
the zer o H amiltonian is in the interior of Ch, because H can simulate any sH' for small s without ancillas (see 
Section III). Thus V Dh' y^ 0, the optimal solution is a boundary point of Ch- The problem of efficient or optimal 
simulation can be rephrased: 

Given H, let Ch be the diagonal subset of the convex hull of {RDhS : R, S £ S0(3)}. Then H' can be 
efficiently simulated by H if and only if Dh' £ Ch- For any H' , sh'\hDh' , which represents the optimal 
simulation, is the unique intersection of the semi-line XDh' (A > 0) with the boundary of C^. The optimal 
protocol can be obtained by decomposing sh'\h Dh' in terms of the extreme points of Ch- 

Since each point in Ch can be decomposed as a convex combination of the extreme points of Ch, each efficiently 
simulated Hamiltonian can be identified with a simulation protocol and vice versa. We will refer to elements in Ch 
as Hamiltonians or simulation protocols interconvertibly 



Central to our problem is the structure oiCn- We will show in Sec. VD that it is a simple polyhedron that we call 



Vh- Its set of vertices, P24, is a subset oiCn containing 24 elements. They are obtained from Dh by permuting the 
diagonal elements and putting an even number of — signs. More explicitly, these elements are ■KiDn'T^iSj, where 



T^o = 1 1 ""i = 



So 



si 



-1 0' 






001' 






"0 


1 1 






"0 1 o' 




"0 r 


1 


, TTa = 


0-10 


, 71'3 = 


1 


, TT4, = 


1 


, TTs = 


1 


1 




1 




0-1 




1 




1 


1 


- 


"-10 ' 




"-1 0" 




0-10 


, S2 = 


1 


, .53 = 


0-10 




0-1 




0-1 




1 





The transformation Dh — > TTiDHi^lsj is physically achieved by iJ — > {Ul.® UsjUjri) H (Ut,.® UJUlJ, where f/^^ = 

(aj + <7k)/V^ for i — 1,2,3 and i,j,k distinct, 1/,^^ = cosjir/S)! ± isin(27r/3)((Ta; + ay + az)/V^ for i = 4, 5, and 
Usi ~ Oi for i — 1,2, 3. These can be verified using Eq. (pi|). The fact that P24 is the set of extreme points of Ch 
means that any optimal simulation protocol only involves the transformations Dh — > T^iDnT^iSj- 

In the next few subsections, we investigate the geometry of Vh, prove that Ch — Vh, and find the optimal solution 
for any H' using the fact Ch — Vh- Then, we restate the solution in terms of a majorization-like relation. 



C. The Polyhedron Vh 



Since P24 and Vh consist of diagonal matrices only, their elements can be represented by real 3-dimcnsional vectors. 
The defining characterization of Vh is the polyhedron with 24 (not necessarily distinct) vertices that are elements of 
P24- We now turn to a useful characterization oiVn as the region enclosed by its faces, 



(x,j/,z)eT'Hiff ■ 



\x\ <hi ,\y\ <hi ,\z\ <hi 
-(I-2/13) <+x + y + z<l 
-(1-2/13) <-x~y + z<l 
-(1 - 2/13) <+x-y-z<l 
-(1 -2/13) < -x + y- z<l 



(27) 



where the fact that H is in normal form, ^13 > 0, and that hi + h2 + h^ = I are used to replace the bounds ^^ hi and 
— (^j hi — 2 mini hi) by 1 and —(1 — 2/13) in Eq. (P7|). Equation (Ej) can be used to determine whether a point, as 
specified by its coordinates, is in Vh or not. The validity of Eq. (|2^) can be proved by plotting P24 (and therefore 
Vh) and verifying that the faces are as given in Eq. ([27|). We first plot Vh for the simple case (hi, /12, /13) = (1, 0, 0), 
for which P24 has 6 distinct points: (±1,0,0), (0,±1,0), (0,0, ±1) and Eq. {^ holds trivially: 




Now, we plot Vh for the most complicated case, hi > h2 > h^ > in Fig. (^): 



(28) 



z = hi 



X — y -\- z = 1 — 2\h_ 




x + y + z = l -2|/i3| 



x + y-z = l-2\h^ 



The equations for 
the faces in the 
background are given 
in boxes. The empty 
faces are given by 
double arrows. 



(29) 



As in Fig. (Eq), Fig. (E9|) is viewed from the direction (1, 1, 1). Three faces are removed to show the structure in the 
back. There are 3 types of faces. There are 6 identical rectangular purple faces on the planes x — ±hi, y — ±hi, z = 
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±/ii. There are two groups each consists of 4 identical hexagonal faces. The first group of 4 consists of the 3 light blue 
faces in the back, and the light blue face in the front. These are the truncated faces of the original octahedron, lying 
on the planes x + y + z — 1, ~x + y — z = 1, —x — y + z = l,x — y — z = 1. The second group consists of the 3 empty 
faces in the front, and the white face in the back. They are inside the original octahedron and are parallel to the 
original faces. They lie on the planes —x — y — z = 1 — 2/13, —x + y-\-z — 1 — 2/13, x — y + z = 1 — 2/13, x-\-y — z= 1 — 2/13. 
Note that each hexagon in one group has a parallel counterpart in the other group. All together, there are 7 pairs of 
parallel faces, each pair bounds one expression in Eq. (p7|). It is straightforward to verify the diagram and Eq. (|27|). 

The plots for other cases, such as when /13 = or /ii = /12, can be likewise obtai ned and Eq. ( p7[ ) be verified. These 
are generally simpler then Fig. (^^, and may admit simpler solutions in Section VE. However, we leave the details 
to the interested readers and move on to prove that Ch = 'Ph- 



D. Proof of Ch ='Ph 



We now show that Ch ~ Vh- Recall that Ch consists of Hamiltonians that can be expressed as Dh' = 
^iPi RiDnSf (by putting s = 1 in Eq. ( p5| ) and using Sj in place of Si). The fact that Dh' is diagonal im- 
plies that only the diagonal elements in each RiDnSf contribute to Dh'\ it is possible for an individual RiDnSf to 
be off-diagonal, but the off-diagonal elements have to cancel out in the sum. To show that Ch = Vh, it suffices to 
show that the diagonal part of each RiDnSf is in Vh, because any Dh' G Ch will then be in Vh- 

Let us consider the diagonal part of RDhS^, represented as a 3-dimensional vector (31,52,53)- Since Dh = 
diag(/ii,/i2,/i3), 



gi — {RDhS )ii — 2_^ RikhkSf.^ — 2_^ RikSikhk 



(30) 



The vectors (hi, /i2, h^) and (51,52,53) are linearly related by 



51 




'hi' 


52 


= R.*S 


h2 


53 




hs 



(31) 



where .* denotes the entry- wise multiplication of two matrices, also known as the Schur product or the Hadamard 
product. It is useful to expand g,; in Eq. (pQ) explicitly 



5i — ^il *5'il hi + Ri2 Si2 /I2 + Ri3 Si3 /13 . 

Then, we can prove the first group of inequalities 

\gi\ < \RtiSti\ hi + \Ri2Si2\ /i2 + |i?i3'5'i3| ^3 < max/ii = hi 



(32) 



(33) 



We use the fact that R,Sg S0(3) to prove the second inequality: R,S consists of orthonormal rows and columns. 
Hence, {\Rii\, \Ri2\, iRisl) and {\Sii\, \Si2\, \Si3\) are unit vectors, and their inner product \RiiSii\ + \Ri2Si2\ + \Ri3Si3\ < 
1. We refer to this argument, which is frequently used, as the "inner product argument". The second group of 
inequalities can be proved by 



E l5«l = E I E ^'fc^'fc ^fc I ^ E ( E \R^k\\S^k\ ) \hk\ <Y.hk^l 



(34) 



The second inequality in Eq. (p3) is due to ^,- |-Rifc||S'ifc| < 1, obtained again by the inner product argument. This 
proves all of 



51 + 52 + 53 < 1 , 51 - 52 - 53 < 1 , - 51 + 52 - 53 < 1 , - 51 - 52 + 53 < 1 



Finally, 



5i + 52 + 53 



i?ll ^ii 
-|-i?21 5'21 I hi 

-R31 S31 




R12 S12 \ I i?l3 5'i3 

R22 S22 \h2+\ +R23 S23 I /13 = ^1^1 + ^2/12 + hhs 

R32 S32 J \+^33 »S'33 ^ 



(35) 



(36) 
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where Xi is the coefBcient of hi in the parenthesis. The inner product argument iniphes |Ai| < 1. Moreover, we will 
prove ^^K > —1 shortly, which implies 

51 + .92 +93> Al/ll + A2/l2 + (-1 - Ai - A2)/l3 

= Ai(/ii - hs) + \2{h2 - hs) - hz 

> ~hi - h2 + hs = -{1 - 2h3) (37) 

where Eq. ( pTI) is the minimum of the previous line, attained at Ai = A2 = — 1 and A3 = 1. We now prove ^^ Xi > —1. 
First, 

/ ^ A,; = i?ii5'ii + i?21'S'21 + R31S31 + R12S12 + R22S22 + R32S32 + ^IS'S'ia + -R23'S'23 + -R33'S'33 = tr(i? S) (38) 

i 

As R, S E S0(3), R^S G S0(3). Each S0(3) matrix is a spatial rotation, therefore having the eigenvalue +1 that 
corresponds to the vector defining the rotation axis. Moreover, any SO (3) matrix has determinant 1. Therefore, the 
eigenvalues are generally given by 1, e**^, e~*"^ and the trace is 1 + 2 cos(0) > — 1. This completes the proof of Eq. ([3^). 
The last 3 of the 4 inequalities 

+ ,91 + .92 + .93 > -(1 - 2/13) , + .91 - 52 - .93 > -(1 - 2/13) , 

-.91 +.92 - 53 > -(1-2/13), - 51 -52+. 93 > -(1-2/13), (39) 

can be proved similarly. For example, consider 

^11 'S'll \ / ^12 512 \ / ^13 '5'i3 \ 

51 - 52 - .93 = ( --R2I 5'21 pll + -^22 5'22 I /l2 + --R23 5*23 /13 ■ (40) 

-R3I 5*31 J \—R32 S32 J V^^33 '5'33 / 

The previous argument for gi + 92 + 53 applies by redefining i?, to be 0—1 x R. 

Altogether, the inequalities in Eqs. (|3Jj), (p4|), and (|39| ) satisfied by (51,52,53) are precisely the defining inequalities 
for Vh in Eq. (p7|). Therefore, the diagonal part of any RDhS^ is in Vh, and Ch = Vh- 

E. Optimization over Vh 

Having proved Ch = Vh, we can solve the optimal simulation proble m gi ven Dh and Dh' by finding the unique 



1 





0" 





-1 











-1 



intersection of the semi-line XDh' with the boundary oiVn (see Section |V B| ) . We now explicitly work out sh\h', i-G- 
the value of A in the intersection, as a function of H and H' . 

Let all the symbols be as previously defined. The intersection is given by u = SH'\H{h'i, h'2, /13), so that 

INIi _ INIi ,,,. 

where \\v\\i for a vector v is the sum of the absolute values of the entries. The set Vh has only 3 types of boundary 
faces. Therefore, there are only 3 possibilities where the intersection can occur: 

1. On the group of faces given by a; + j/ + z = 1, —x + y — z — 1, — x — y + z — l,x ~ y — z = 1. In this case. 



V 



|i = 1, and sh'\h 



1 



2. On the group of faces x + y — z = l — 2/13, x — y + z = l — 2/13, — x + y + z = 1 — 2/13, — x — y — z = 1 — 2/13. In 
this case, ||tr||i = 1 - 2/13, and sh'\h = h{+hl+\h'J - 

3. On the group of faces x = ±hi,y = ±/ii,z = ±/ii. In this case, v = ^{h'i,h'2,h'^) (note h[/hi > 0), 
||i7||i = ^{h'l + /i2 + I/13I) (not constant on the face), and sh'\h — ^■ 

Note that when H' is in normal form, v can only fall 0'iix + y + z = l, x + y — z— 1 — 2/13, and x — hi in each of 
case 1, 2, and 3. We now characterize the {h[,h2, /I3) belonging to each case. 
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Case 1. Note that the face of Vh ona; + i/ + z = lis the convex hull of (ft.i,/i2,ft-3) and all permutations 
of the entries. The hexagon contains exactly all vectors v majorized by (/ii, /12, ^3), v -< (/ii, /12, /is) (see next 
section for definition of majorization). Hence, {h'i,h'2,h'^) is in case 1 if and only if it is proportional to some 

V'< (/ll,/l2,/l3)J] 



• Case 3. In this case, v 
with vertices (/i2, /is), (^3, ^^2), (-/12, -h^), (-/13, -/i2) 



/ 7 All hr, h-[ Alo 



Thus {h'l, ft,2, /13) is in case 3 iff (-^, -^7^) is within the rectangle 





(/!3,ft2) 


A 


^^(/!2,/»3) 


(-fts.-fta)^^ 


|j/ + z| < /j2+/l3 


(-fts-ftsf 


|?/-2| < /J2-/J3 



(42) 



Hence, (ft'^, /I'j, /13) is of case 3 
•ff ^1^2 

iff 



K 

hi 



< /i2 + /13 and 



< 



K + K 



and 



hlhn ft-l/lQ 



hi 



h[ 



<h2 



h?, - h-^ 



< 



h'l 



h'^ - h'. 



(43) 
(44) 



• Case 2. This contains all {h'l, h'2,h'^) not in case 1 or 3. 

The intersection on a boundary face can be easily decomposed as a convex combination of at most 3 vertices in ^24- 
The decomposition directly translates to an optimal protocol (using the discussion at the end of Section V_B) with at 
most 3 types of conjugation. 



F. Optimal simulation, polyhedron Vh, and s-majorization 

The problem of Hamiltonian simulation can also be analyzed from the perspective of a majorization-like relation, 
which provides a compact language to present the main results of this paper. 

Let us recall the standard notions of majorization and submajorization as defined in the space of n-dimensional 
real vectors. Let u be an n-dimensional vector with real components Ui, i = 1, • • • ,71. We denote by u^ the vector 
with components wj > Uj > • • • > u^, corresponding to \ui\ decreasingly ordered. Then, for two vectors u and v, u is 
submajorized or weakly majorized by v, written u -<w v, if 






u{+U2<v[ + V^ 



2 1 



+ M2 H u]^<vi+v'2 + 



(45) 
(46) 

(47) 
(48) 



In case of equality in the last equation, we say that u is majorized by w, and write u < v. 

These notions can be extended to real matrices. Suppose M and N are two nxn real matrices. Let sing(M) denote 
the set of singular values of the matrix M (and similarly for TV). Then, M -<yj N, when sing(M) ^^i sing(iV). Thus, 
majorization endows the set of real matrices with a partial order, and a notion of equivalence. 



M - O1MO2 



(49) 



The fact (/li, h'2,h'^) -< {hi, /12, /ig) is a necessary condition for efficient simulation is independently proved in Ref. 
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where Oi G 0{n) are orthogonal, and the transformation M -^ O1MO2 preserves the singular values. A "convex 
sum" characterization of weak majorization 

M ^^ N ^ M = Y,P^O^lNOa (50) 

i 

also holds |g6[, meaning that N always weakly majorizes a (left and right) orthogonal mixing of itself. 

We now introduce the notion of special majorization, s-majorization for short, on n x n real matrices M . First, we 
define a new equivalence relation in terms of any R,Sg SO(n): 

M -^ RMS . (51) 

Then, we introduce the s-majorization relation by means of the "convex sum" characterization, 

M ^s N -^ M = ^piRiNS^ . (52) 

i 

That is, M is s-majorized by N when M is a (left and right) special orthogonal mixing of TV. Equations (|5|) and 
(b2) suggest s-ordering for a real vector u, u^'^: 

{\u\[Mi,---Ml_^,sg{U,u.)\u\i). (53) 

In other words, we rearrange the absolute values of Ui in decreasing order, and put the sign of the last element to be 
the product of all the original signs. Our results in Section VH-VE imply, for n — "i, the following characterization 
of s-majorization: 

Let u and v be 3-dimensional vectors, and their s-ordered versions be u-^' — {ui,U2,U3) and w^"' = 
(vi,V2tV3). Then u ^g v ii and only if 

Ul<Vi, 
Ml + U2 - U3 < Vl + V2 - vs , 
Ui + U2 + U3 < Vi + V2 + V3 . (54) 

Remark Note that u ^s v ^ u -<w v. Moreover, when sg{IliUi) = sg{JIiVi) and ^^ Ui = ^^ Vi, then -K^, -<s, and -< 
are all equivalent. 

We can now state our result in Hamiltonian simulation in the language of s-majorization: 

Theorem: Let H ~ J^i ^t'^i ® ^i ^'^d H' ~ ^. h[ai (g) ai, h — {hi, ft,2, ft.3), and h' — {h[, h'2, /13). Then 

H' <Lu H ^ h' <, h. (55) 

The optimal simulation factor is given by Sfjnu = max^/i'^^h s. 

Proof: Note that when both h and h' are s-ordered, Eq. (|2j) that characterizes Vh reduces to 

K < hi , 
h[ +h'2-h'^<hi+ /i2 - hs , 
/I'l + h2 + h'3<hi+h2 + /i3 . (56) 

which is the condition for s-majorization. Thus h' E Vh iff h' ^s h, and H' <lu H <^ h' ^s h. O 

VI. HAMILTONIAN SIMULATION WITH LU+ANC 

In this section we will show that the use of uncorrelated ancillas does not help when simulating one two-qubit 
Hamiltonian with another, so that all results on efficient and optimal simulation under LU hold under LU-|-anc. We 
prove this by describing the most general LU-|-anc protocol and reducing it to an LU protocol. 

In this scenario, qubits A and B are respectively appended with ancillas A' and B' , which have finite but arbitrary 
dimensions. The initial state of A'B' can be chosen to be a pure product state |0^') \0b')- At the final stage of 
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the simulation, the ancillas A' and B' may be correlated, but A'B' is uncorrelated with AB if the latter is to evolve 
unitarily according to H'. The local unitary transformations Ui and Vi can act on A A' and BB' respectively. This 
feature distinguishes LU+anc from LU. 

The most general LU+anc protocol to simulate H' with H can be described as 



(^Uf<»VfxUn<»Vn e-'"'" Ul^V^x---xUi 



Vi e 



'"'' Ul (g) V^ ) IV') «) \0a') «) |0i3' 



{e-'" ' \^) ) ® {WA'B'ihr ■ ■ ,t„)\OA') ®\0b') ) ■ (57) 



In Appendix O we have shown that for infinitesimal times Eq. (p7|) leads to 



sH' 



AB = {Oa'\® {Ob'\ [ J2pk Uk ®VkiH® Ia'B') ul ® Vl ] |0a') ® |0 



where pk = tk/t and s = t' /t. Let Mk = {^A'Wk and N^ = (0_B'|Vfe. We can write Eq. ( p8| ) as 

sH' = J2pk Mk ^Nk{H® I A'B') Ml nI . 



(58) 

(59) 
I A'B' and the local 



Note that this is the LU+anc analogue of Eq. (|24| ) for LU. In this case, H is replaced by H < 
unitaries are replaced with more general transformations. 

We focus on just one term in the convex combination of Eq. (|5|), M i^ N (BAB'S) Ia' B') M^ ig) N\ with M = {OA'\Uk 
and N = (Os'|Vfe. We will show how to obtain the same contribution to H' using only local unitaries on A and B to 
establish the equivalence of LU and LU+anc. First, note that 



M(g)N [Hab ® Ia'B')M^ ®N'^ =£ao £b{H) 



(60) 



where £a{t) = M{t S Ia')M^ and similarly for £b- We emphasize that £a,b are linear operators on matrices that 
are not necessarily quantum operations p5[ , despite various resemblances to the latter. One can check that £a is 
unital, i.e. £AiI) = I, by using M = (Oa^. Furthermore, £a is completely positive |g^, because an operator-sum 
representation £a{t) = ^iFixF^ can be obtained by expanding Ia' in terms of some basis {|m')}, and by writing 
Fi = M\iA') = {OA'\U\iA')- However, in general, £a is neither trace nonincreasing or trace nondecreasing, though 
tTAj^i^i Fi = tTAj2ii^A'\U'^\0A'){0A'\U\iA') = 2. For each Fi, we can obtain the singular value decomposition 
J^i = W2iQiWu, where Wu and W2i are unitary, and 



Q^ = 



Qii 

q^2 



is diagonal and positive semidefinite. Altogether, 



£a{t) = Y, liQa + 4) W2^ Q^ Wu r Wj, Q, Wl 



4 



(61) 



(62) 



cos Oi 
sin 6* 
Q,(r) = (1 -cos 61 sin 6*) ItI 



where Qi = \/2 



and cos^i = In/ ^/qfi + Qa- We now claim that we can replace the action of Qi by 
cos sin UzTaz, i.e., replacing £a by the following: 



^A{r)=Y.\{ql 



qh)W2.Q,iWuTWl)W^^ 



(63) 



It is straightforward to verify that 

QiIQi = I + cos20 a^, Qia^Qi 



sin26'crj;, Q iffyQ I = sin 26 ay , Qia^Qi = cos29 I + a^ 



Qiil) ^ I , Qi (cTx )= sin 26* CTa;, Qi{ay) = sin29 ay , Q^{az) ^ a^ 



(64) 
(65) 



As H is purely nonlocal, the input to £a in Eq. (pfl) is always traceless. Now, consider how £a and £a act on a^.y^z- 
The only difference is that Qi is not producing the extra / component produced by QiazQi- This extra / component 
has a final contribution as local terms in Eq. (|6y), which can be ignored. Finally, we note that ^^ ^{ifi + ^iz) — 
^'^itrFj Fi — 1, so that £a is indeed a convex combination of the individual terms, each in turns a mixture of 
unitary operations on A. Applying the same argument to £b, Alice and Bob only need to perform local unitaries in 
the simulation step of Eq. (pO) . 
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VII. DISCUSSION 



First, we point out that the normal form for Hamiltonians acting on 2 qubits (See. V A) is symmetric with respeet 
to exchanging the systems A and B. More formally, let S{Mi (g) M2) = M2 ® Mi be the (nonlocal) SWAP operation. 
Then H =lu S{H). This has important consequence - any task generated by the Hamiltonian can be done equally 
well with the role of Alice and Bob interchanged. 



In higher dimensions, the 
the Hamiltonian (see Ref. |^ 



property H =lu S{H) no longer holds. For example, H ^lu S{H) and S{H) ^lu H for 
J for proof): 



H = 



1 0" 




0-10 


® 








1 
1 
0-2 



(66) 



In fact, ii H = Hi (E) H2 where Hi and H2 are members of a traceless orthogonal basis with different eigenvalues, 
S{H) -^LU H and H ^ S{H). This also has important consequences - in higher dimensions, the nonlocal degrees of 
freedom of a Hamiltonian cannot be characterized by quantities that are symmetric with respect to A and B, such 
as eigenvalues of H (independently reported in Ref. pl[|). Any normal form necessarily contains terms of the form 
Cijrji (g) rjj for some nonzero Cij and the matrix with entries c^ cannot be symmetric. 

Second, we revisit the notion of efficiency in Hamiltonian simulation. Our definition of H' < H depends on the 
normalization of both H and H' . One method to remove the normalization dependence is to require h'-^-\-h'2 + \h'^\ = 1. 
Alternatively, we can consider the product SfJ\H'SH'\H^ that measures the inefficiency of interconverting H and H' 
independent of the normalization of the Hamiltonians. We found that (proof omitted) when /13 > 0, sh\h'Sh'\h ^ s- 
Otherwise, sh\h'Sh'\h > 5, with the lower bound attained at /i = (1/3, 1/3, 1/3) and h' = (1/3, 1/3, —1/3). 

Third, we have considered the optimal simulation of one two-qubit Hamiltonian using another, both arbitrary but 
known. We can apply the characterization of Vh to analyze other interesting problems. For example, inverting a 
known Hamiltonian is equi valent to setting H' = —H. Without loss, assume ft.3 > and hi + h2 + h^ — 1. Using 



the analysis in Section VE, the intersection is of case 2. Therefore, s_h|h = ^(1 ^ 2/13). The worst case is inverting 



[ux ® <Tx ~\- <7y ® Uy + (Jz ® f^z) in which case s_if|_ff = 1/3. In contrast, any protocol for inverting an unknown 



Hamiltonian can invert the worst known Hamiltonian, thus s_ 
protocol: 



-H\H < 1/3. This is achievable using the following 



./ 



-iHt' 



iHt' 



{ax(Jy®I) e •'"* {ayaz®I) 






I 



-i(~H)t'/3 



(67) 



We can also improve on the time requirement for simulation protocols for n-qubit pairwise coupling Hamiltonians m 
with our construction. Instead of selecting a term by term simulation using a single nonlocal Pauli operator acting on 
a pair of qubit, one can directly simulate the desired coupling between the pair with any given one in a time optimal 
manner. 

VIII. CONCLUSION 

We have discussed various notions of Hamiltonian simulation. Focusing on dynamics simulation, we show its 
equivalence to infinitesimal simulation, and the intrinsic time independence of the protocols. We also show the 
possibility of simulating one nonlocal Hamiltonian with another without ancillas in any two d-dimensional system. 
Our main results are on two-qubit Hamiltonians, in which case, for any Hamiltonian H, we characterize all H' 
that can be simulated efficiently, and obtain the optimal simulation factor and protocol. We obtain our results by 
considering a simple polyhedron that is related to some majorization-like relations. Our results show that the two- 
qubit Hamiltonians are endowed with a partial order, in close analogy to the partial ordering of bipartite pure states 
under local operations and classical communication pTJ ]. 

We have restricted our attention to simulation protocols that are infinitesimal, one-shot, deterministic, and without 
the use of entangled ancillas and classical communication. We also restricted our attention to bipartite systems. 
Extensions to the unexplored regime, and alternative direction such as other nonlocal tasks will prove useful, and are 
being actively pursued. 



16 



IX. ACKNOWLEDGEMENTS 

We thank Wolfgang Diir, John Smolin, and Barbara Terhal for suggestions and discussions that were central to this 
work, as well as Herbert Bernstein, Ben Recht, and Aram Harrow for additional helpful discussions and comments. 
CHB and DWL are supported in part by the NSA and ARDA under the US Army Research Office, grant DAAG55- 
98-C-0041. GV is supported by the European Community project EQUIP (contract IST-1999-I1053), and contract 
HPMF-CT-1999-00200. 

APPENDIX A: NOTIONS OF SIMULATION 

We consider various notions of using a Hamiltonian H to simulate the evolution due to H' for time t' . 

In dynamics simulation, the evolution of the system is close to e~*^ * after an operation time of fit" for constant 
H and V t" € [0,i']. It is possible to relax this requirement, so that, /i(t") is a function of t", and without loss of 
generality, /i(t") is nondecreasing. We call this "variable rate dynamics simulation". Finally, in gate simulation, the 
only requirement is that, the final evolution is given by e~*^ * . 

As an analogy, let H' be driving along a particular highway from my house to your house at 100 km/hr. Dynamics 
simulation is like driving, biking or walking along the same highway at any constant speed. Variable rate dynamics 
simulation is like driving along the highway at variable speed, for example, when there is stop-and-go traffic. The 
vehicle is always on the trajectory defined by H' . Finally, gate simulation is like going from my house to your house 
by any means, for example using local roads, or flying a helicopter. 

It is important to note the difference between dynamics simulation (or infinitesimal simulation) and variable rate 
dynamics simulation. For example, iterating infinitesimal simulations to perform dynamics simulation, the ancillas 
are implicitly discarded after each iteration, and new ones be used next. However, it is possible in variable dynamics 
simulation that used ancillas can subsequently be used to accelerate the simulation. Such phenomena are known in 
entanglement generation 0]. The more complicated analysis for variable dynamics simulation will be addressed in 
future work. 



APPENDIX B: SIMULATING ZERO HAMILTONIAN IN D x D WITHOUT ANCILLAS 



In Ref. ||23|] , it is shown that for any d-dimensional square matrix Af, 



(Bl) 
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u^,= 



w'^-i 
and a; is a primitive rf-th root of unity. H can simulate using the protocol: 



which is local and can be removed 



(B2) 



ilKBt 



(B3) 



APPENDIX C: ARBITRARY HAMILTONIAN SIMULATION IN D x D WITHOUT ANCILLAS 

Let H and H' act on two d-dimensional systems. We use the following (nonorthonormal) basis for traceless hermitian 
operators acting on a d-dimensional system: 
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Let H — J2ij Cij ili ^ Vj and H' — ^. ■ c'^jiji ®rij. To show that sH' <lu H for some s > 0, it suffices to show that 
s rji® rji <Lu H, since rji ® 771 =lu i^i ® Vj foi' ^U hj- Furthermore, rji (g) 771 can be simulated if one can simulate 

N)(«l «) \j}{j\ and -\i'){i'\ ® I/) (/I for any i,j,i\j'. 



Without loss of generality, cn 7^ 0. We first use H to simulate its diagonal components, Hd = X) 
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Hd can further be used to simulate cii|2)(2| (g) |2)(2|, using the protocol 
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(C2) 



This corresponds to Alice and Bob each applies an averaging over all the computation basis states except for 1 2) . Since 
all r]i:^i are traceless on the subspace spanned by |i 7^ 2), they vanish after the averaging, leaving only a contribution 
by ?7i (g)?7i: 
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(C3) 



which is equivalent to cn |2) (2| (g) |2) (2| up to local terms. It remains to obtain a term with sign opposite to cn. If some 
Cki y^ has a sign opposite to cn, we can simply repeat the same procedure, with Alice applying an averaging over 
all \i y^ k) and Bob applying an averaging over all |j ^ l). If all Cij has the same sign, Alice can apply an averaging 

over all \i ^ 1) and Bob can apply an averaging over all \j ^ 2) to obtain — YliZi Cii|l)(l| ^ |2)(2|, completing the 
proof. 

APPENDIX D: INFINITESIMAL SIMULATION AND TIME INDEPENDENCE 

The most general simulation protocol of H' with H using LU+anc can be described by 

( [// (g F/ X [/„ (g K e'^^*" C/^ (g y„^ • • • C/i (g VI e-^^*i C/| ®V^){ \^) (g \Qa') ® \Qb') ) 



( 



-iH't' 



\iP))®{Wa'b' \Qa')®\Qi 



(Dl) 



where the equality must hold for all possible states |^) of system AB. Here the unitaries Ui and Vi, acting on AA' 
and BB' respectively, and the partition {ij} of the time interval t = J2i^i' correspond to all the degrees of freedom 
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available for the simulation of H' for time t' . The initial states of the ancillas A' and B' is \0a') "8) \0b'), and Wa'B' is 
their residual, unitary evolution, which is determined by the other degrees of freedom and may create entanglement 
between A' and B' . 

We have argued earlier that optimal dynamics simulation can always be achieved by a protocol for simulating 
infinitesimal evolution times f . This also implies t being infinitesimal. Recall that pi = ti/t and s = t' /t. We can 
expand Eq. ( |D1| ) to first order in t to obtain 

Uf®Vfy.\l^ itY^P^ U, ®V,{H® Ia'b') UJ V^^ 1 |0^.) |0b') = ( / - ttsH' ) ® ( Wa'b' \0a') <» \0b') ) ■ (D2) 



The validity of Eq. (Dl) for all |?/;) is used to obtain Eq. (D2), each term of which is taken to be an operator on AB. 
It follows from Eq. (D2) that 



which implies that 



Uf\OA') ) ® ( Vf\OB') ) - /as ® ( Wa'B'\Oa') ® |0i3') ) + 0{t) , 

Uf \0a') =Ia®{ Wa' \0a') ) + 0{t), 

Vf\0B')=lB(^{WB'\0B'))+Oit), 

Wa'B' \0a') ® \0b') = ( Wa' \0a') )®{Wb' \0b') ) + 0(t) ■ 



(D3) 

(D4) 
(D5) 
(D6) 



Equation (D6) implies Wa'B' is a product operator to zeroth order in t. Redefining Uf and Vf is n e cess ary, we can 



assume Wa' = Ia' and Wb' = Ib' ■ Explicitly writing down the most general 0{t) terms in Eqs. (D4)-(D£), we obtain 

Uf \0a') = ( Iaa' ~ itKAA' ) |0a') + O(t') , (D7) 

Vf \Qb') = ( Ibb' - HKbb' ) |0b') + 0(t2), , (D8) 

Wa'B' |0a') ® \0b') = ( [Ia' - itKA')\OA') ) ® ( {Ib' - itKB')\OB') ) - UKa'B' \0a') ® |0b') + ©(i^). (D9) 

where the unitarit y o f the operator s on the LHS implies the hermiticity of Kaa', Kbb', Ka', Kb', and Ka'B'- 
Substituting Eqs. (^-(^) in Eq. (Q imphes 



sH'®{Ia'B'\0a')®\0b')) = 

{J2p^U,®V,iH®lA'B')UJ 



y/ + Kaa' + Kbb' - Ka' - Kb' - Ka'b' ) |0^ 



|0i 



0{t). 



(DIO) 



Projecting this equation on the left onto {Oa'\ ^ {Ob'\, the terms Kaa' + Kbb' — Ka' — Kb' — Ka'B' become local or 
identity terms. Taking into account that H' has zero trace and no local terms (recall Section I.C.I), their contributions 
vanish, and we obtain 



SH' = (Oa'I ® (Ob'I [Y,P^U^®V^{H® Ia'B') U} ® v} ] |0a') ® |0s 



In the case we do not have ancillary systems, Ui and Vi only acts on A and i3, and Eq. (Dll) reads 

sH' ^^p,U^®V^HU\® V} , 



(Dll) 



(D12) 



In Eqs. ( |D11| ) and ( p312| ), the dependence of the equation on the original infinitesimal times t and t' is only through 
s ~ t' /t. This implies any protocol for t and t' , applies to at and at' within the infinitesimal regime. Thus the 
protocol, namely, the set {Ut, Vi,pi} can be considered being independent of t in the infinitesimal regime. 
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We show that S{H) ^lu H, which also implies H ^lu S{H). We can renormalize both H and S{H) so that the basis is 

orthonormal. The condition S{H) <lu H implies 

H2®Hi= ^p,{U,mul) ® {ViH2V^) 

for some unitary Ui,Vi. Multiplying the product of the above with H2 ® H\ and taking the trace, 

tr{Hl)tr{Hl) = ^pMH2UzHiUl)tr{HiV^H2V^) 

i 

However, tr{H2UiHiUl) < tr(_H'|) and tr{HiViH2V^^) < tr{Hl), with equalities hold only when UiHiU} = H2 and 
ViH2V^ = Hi, implying Hi and H2 need to have the same eigenvalues. It is well known that, in d-dimensions, the traceless 
matrices diag(l, 1, ■ ■ ■ , 1, — fc, 0, 0, • ■ ■ , 0) are part of a traceless orthogonal basis for hermitian d x d matrices (for example, 
the Gell-Mann matrices) completing the proof. 
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